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Abstract 

Quantum nonrelativistic systems with 2x2 matrix potentials are investigated. Phys- 
ically, they simulate charged or neutral fermions with non-trivial dipole momenta, inter- 
acting with an external electric field. Assuming rotationally invariance of the Hamiltonian 
all such systems allowing second order integrals of motion are identified. It is shown that 
the integrals of motion can be effectively used to separate variables and to reduce the 
systems to decoupled ordinary differential equations. Solutions for two of the discussed 
problems are presented explicitly. 
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1 Introduction 



Exactly solvable systems of quantum mechanics are favourite subjects for many physicists 
and mathematicians. The beauty of such systems (like the Hydrogen atom or the harmonic 
oscillator) is that they are simple enough to be solved in a way free of uncertainties and 
inconveniences of various approximate approaches. Meanwhile, they are sufficiently complicated 
to model the physical reality. In addition, the complete sets of their exact solutions supply us 
by convenient bases for expansions of solutions of other problems. Many of exactly solvable 
systems admit nice hidden symmetries which are very interesting on their own account. 

The exact solvability of quantum mechanical systems is usually caused by their specific 
property called superintegrability. The system with n degrees of freedom is integrable if it 
admits n — 1 commuting integrals of motion in addition to its Hamiltonian. The system is 
superintegrable if it is integrable and admits at least one more integral of motion. 

The systematic search for superintegrable systems started with paper p. We will not 
recount the details of rather inspiring history of this search and discuss all obtained fundamental 
results related to quantum mechanical and classical scalar systems. In contrary, we restrict 
ourselves to discussion of superintegrable systems with spin. 

Such systems where studied methodically in recent papers The problem of classifi- 

cation of superintegrable systems with spin was stated in [2] where 2d systems with spin-orbit 
interaction and first order integrals of motion had been presented. Then these results were 
generalized to the cases of three-dimensional Euclidean space [3] and second order integrals of 
motion |1] . The results presented in [1] were restricted to rotationally invariant systems and to 
integrals of motion which are rotational vectors or scalars. 

However, the spin-orbit coupling is not the only spin effect which can be present in quantum 
mechanical systems. One more coupling which is very important and has a well observable effect 
is the dipole, or Pauli interaction. This interaction is represented by the Stern-Gerlach term 
~ S • B (where S and B are the spin and magnetic field strength vectors) ore, more generally, 
by a matrix term linear in an external field. Moreover, the Pauli interaction affects even neutral 
particles provided, like neutron, they have non-trivial dipole moments. 

An important example of a 2d superintegrable system with dipole spin interaction was 
presented long time ago by Pron'ko and Stroganov [5j . However, a systematical search for such 
systems is only at the beginning. The classification of 2d systems with the first order integrals 
of motion was presented in [B] , while the 3d systems with Fock type dynamical symmetry were 
derived in [7]. Generalizations of the Pron'ko-Stroganov system to the case of arbitrary spin 
were discussed in |8] and |9]. However, a classification of integrable and superintegrable systems 
with the dipole interaction, admitting higher order integrals of motion, was not carried out till 
now. 

In the present paper 3d superintegrable systems with dipole interaction are classified. We 
restrict ourselves to systems which are invariant w.r.t. the rotation group and admit second 
order integrals of motion. The list of such systems is not too long, but two of them are defined 
up to arbitrary functions, and so the number of inequivalent systems is infinite. In particular, it 
includes a system with Fock symmetry discussed recently in [7], and a supersymmetric system 
with a matrix potential. For a classification of supersymmetric matrix potentials see papers 

m-m- 

Superintegrable systems are in many cases exactly solvable, see, e.g., [I3]-[I6]. In particular, 
all the spinless 2d superintegrable systems whose integrals are given by second order differential 
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operators are exactly solvable as well, and it was conjectured that this property is valid also 
for higher dimensional superintegrable systems jT5] . 

For systems with spin the situation is more complicated due to the presence of the additional 
dichotomous variable. As it was noticed in [7J, such systems with 3 spatial variables are not 
necessary exactly solvable even if they admit more than 3 second order integrals of motion. 

Thus the solvability of the discussed systems should be examined separately. Of course, it 
is impossible to solve all systems with arbitrary functions, presented in the current paper. We 
restrict ourselves to two particular cases which are seemed to be physically interesting since 
include the field of point charge. The corresponding exact solutions are presented in section 6, 
where some other exactly solvable problems are indicated also. In addition, we apply integrals of 
motion to reduce generic (i.e., including arbitrary functions) eigenvalue problems to decoupled 
systems of ordinary differential equations for radial wave functions. However, these equations 
in general are not exactly solvable. 

2 Schrodinger-Pauli equations for neutral particles 

To describe spin effects in non-relativistic quantum mechanics the Schrodinger equation should 
be generalized by Pauli term proportional to the scalar product of spin with vector of the mag- 
netic field strength. For neutral particles with non-trivial dipole moments (e.g., for neutrons) 
this term becomes dominant since in this case the minimal interaction is absent. The Pauli- 
like term is requested also for description of interaction of charged particles having non-trivial 
electric moments with an external electric field. 

In other words, there are several reasons to study the Schrodinger-Pauli equations of the 
following generic form: 



where cr is the matrix vector whose components are Pauli matrices, K and V are vector and 
scalar external fields. Moreover, A and u are coupling constants and E denotes an eigenvalue 
of Hamiltonian H. 

Since we are interested in systems with a spin dependent interaction, constant A is supposed 
to be nonzero. To obtain more compact formulae in the following calculations let us rescall 
variables and reduce the Hamiltonian to the following simplified form: 



Just Hamiltonians ([2]) which include the generic 2x2 matrix potential V will be the subject 
of our classification. Our goal is to find all possible external fields F = {F^,F'^,F^) and F^ 
such that the systems with such Hamiltonians be superintegrable. 

Notice that symmetries of systems with 2x2 and 3x3 matrix potentials have been inves- 
tigated in papers [17] and [H] . However, it was done only for diagonal potentials V depending 
on time and one spatial variable. 




(1) 




(2) 



where A is the Laplace operator and V^(x) is a matrix potential . For this purpose we change 
in (HD 

E ^ E = 2mE, K ^ F = AK, V ^ F° = 2muV. (3) 
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3 Determining equations 



Let us search for the first and second order integrals of motion for systems described by equation 
([1]). By definition such integrals of motion Q are first and second order differential operators 
commuting with Hamiltonian H: 

[H, Q] = HQ-QH = 0. (4) 

We suppose these operators be formally self adjoint. Then, without loss of generality, they can 
be written in the following form: 

Q = V4, V J + ia'^{A^^ VJ + a'^Q" (5) 

where A'"' and il'" are (unknown) real functions of x, {^^^"■'', Va} = $'^''^Va+Va$^"^, Va = 
are Pauli matrices with a" being the 2x2 unit matrix. In addition, here and in the 
following the summation is imposed over the repeated indices. Moreover, in all equations the 
Latin and Greek indices take the values 1,2,3 and 0,1,2,3 correspondingly. 
Substituting ([2]) and ([5]) into (j4]), using the relations 

a a = + le a 

where e"^'^'^ is the Levi-Civita symbol, and equating coefficients for linearly independent matrices 
and differential operators, we obtain the following system of determining equations for functions 
(^^^ab^ ^i^a^ and F^: 



$f ^ + $f + $f = 0, (6) 

Ar + = 0, (7) 

Ar + + e'^"'^^^*""'' = 0, (8) 

^Oakpb ^ ^bakpO _ 2^bckp^capk + = Q, (9) 

^akpO ^ ^Okpa ^ ^ack^cpk ^ ^nka^ncdpk^ ^ (^q) 

^/.afc^M + f]0 = 0, (11) 

^^,kp^^ = 0. (12) 

Here the subindices denote derivatives w.r.t. the independent variables, i.e., = etc. 



Thus to classify Hamiltonians (|2]) which admit first- and second-order integrals of motion 
we are supposed to solve the system of determining equations (l^- ffT^ . Moreover, searching for 
the first order integrals of motion we should a priori set ^^""^ = 0. 

4 Classification results 

The system fP)- f[T^ is rather complicated and includes 77 coupled nonlinear partial differential 
equations for 44 variables. However, it is rather symmetric, and some of its constituents are 
easy integrable. The general solution of this system with two independent variables for the case 
= was found in [7]. 
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In this paper we find solutions of equations (16|)-( !T2|) compatible with the supposition that 
Hamiltonian is invariant with respect to the rotation group 0(3). In other words we suppose 
that H commutes with generators of this group which are nothing but components of the total 
orbital momentum vector J: 

J = xxp + ^cr. (13) 

This condition reduces the general form of external fields to: 

F° = 0(x), F'^ = xV(x) (14) 



where and are functions of x = a/x^ + + x\. 

The system flBl)- ffT^ with additional conditions f[T^ is algorithmically solvable, see Ap- 
pendix. As a result we obtain the following list of Hamiltonians ([2]) together with the admissible 
constants of motion additional to (fT3l): 



72 ■ ■ X 



Qi = cr ■ L + 1 + Act • n; 



i7 = i72 = -V2 + ^-nr + f --, 



Q2 = (io- ■ p + /) (o- ■ L + 1) + -cr ■ n; 



(16) 
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H = H^ = -V^ + -cr ■ n + 



X x^ x' (17) 

Qiy Q3 = Q2I /=a; 



^ 1 , , , , Axcr ■ X ^ ^ 

R = -ipy< J - 3 y<p) + —^^ — . 

Here L = x x p is the orbital momentum, n = -, <f{x) and / are arbitrary function of x, 
f = ^ 

Hamiltonian f lT7|) is a particular case of operators presented in f lTSjl and f[T6|) . which corre- 
sponds to a more extended number of integrals of motion. 

Hamiltonian f lTS]) admits three integrals of motion, additional to f[T^ . They are components 
of vector R, which generalizes the Laplace- Runge-Lenz vector to the case of a system with spin. 

The system with Hamiltonian (|T8|) possesses the dynamical symmetry w.r.t. group 0(4) 
whose generators are vectors J and R. This system was discussed in paper \^ where its shape 
invariance was proven and exact solutions of the corresponding eigenvalue problem ([1]) where 
found. 

Notice that all Hamiltonians (IT^ - (IT5]) admit one more (discrete) symmetry. Namely, they 
commute with the following operator 

Qi = {(T-L + l)p (19) 
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where p is the space inversion operator which changes the sign of independent variables, i.e., 
acts on the wave function as pipi'x.) = ■?/'(— x). 

One more symmetry which includes the space reflection is valid for Hamiltonians H2. 
Namely, these Hamiltonians commute with the following operator: 

Q5 = {a- ■ p + f)p + acr ■ n. 

Of course, this symmetry with / = - is valid for Hamiltonian H^. 



5 Algebraic properties of integrals of motion 

Let us discuss some general properties of the constants of motion presented in the previous 
section. 

By construction, operators Qi, Q2, Q3 and R commute with the corresponding Hamilto- 
nians, i.e., satisfy conditions (jl]). 

Operators Qi, Q2, and Q3 are rotational scalars and so commute with the total orbital 
momentum (fT3l): 



[Q„J] = 0, a = 1,2, 3. (20) 
In addition, these operators satisfy the following algebraic relations: 

Ql = 3' + \' + \, (21) 

QI=(3' + ]]h2 + ^, (22) 
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^2 / t2 



1\ „ «2 



Qi = l^J' + l)H3 + (23) 

QiQa + Q3Q1 = «A, [Ql,Q,] = 0, [QlQi] = 0. (24) 

We see that the integrals of motion form rather non-trivial superalgebraic structures. Re- 
lations (123!) and fl24l) will be used to explain the degeneration of spectrum of Hamiltonian 

In the particular case ip = integrals of motion f ll3p . f llSp together with operators D = 
Xipi + X2P2 + X3P3 and K = form a basis of the seven-dimensional Lie algebra since the 
following commutation relations are satisfied: 

[H,,D] = -2iH, [K,H,] = iD, 

[r, /] = ie'^^V^ [K, D] = 2iK ^ ' 

while all the other commutators are trivial. In other words, there is a direct sum of the conformal 
algebra so(l,2)9 {Hi, D, K), the Lie algebra of the rotation group so(3) 3 { J^, J^, J^) and 
the one-dimension algebra spanned on Qi. Thanks to the conformal symmetry the discussed 
system can be interpreted as a model of conformal quantum mechanics, see [19] for definitions. 

Finally, consider the superintegrable system fllSp admitting vector integrals of motion. Com- 
ponents Ra of vector operator R satisfy the following commutation relations: 

[Jay Jb] '^^abcJci \Ray Jb\ '^^abcRci (26) 
[Ra, Rb] = —^lEabcJcHi- 
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Being considered on eigenvectors of Hamiltonian corresponding to coupled states, algebra 
(12^ is isomorphic to the Lie algebra of group 0(4). In other words, the system with Hamiltonian 
H4 admits the same dynamical symmetry as the Hydrogen atom. The detailed analysis of this 
system is presented in [7]. 



6 Exact solutions 

Since Hamiltonians Hi and H2 are defined up to arbitrary functions, they represent an infinite 
set of superintegrable models. In this section we consider two important particular cases of 
such models which involve the Coulomb potential, and present a constructive way for finding 
solutions for the systems with arbitrary potentials. 



6.1 Charged particle with electric dipole moment interacting with 
the field of point charge 

Let us start with Hamiltonian Hi. It includes an arbitrary scalar potential and a dipole 
interaction term ~ cr ■ F with coupling constant A and external field F = p-. This field can 
be interpreted as an electric field generated by a point charge. Thus it is naturally to choose 
ip = ^, then the corresponding operator Hi can be interpreted as a Hamiltonian of a charged 
particle with spin 1/2 and a non-trivial dipole electric moment. 

Consider the eigenvalue problem for such specified Hamiltonian Hi. 

f-A + --)^ = Ei,. (27) 

Here = ipi'^) is a two-component function which is supposed to be normalizable and vanishing 
at X = 0. In addition, to obtain a system with coupled states we suppose that a > 0. 

Equation (^7\i admits three constants of motion J3, and Qi which commute each other. 
Thus we can expand = V'(x) via eigenvectors Qj^i^^^{ip, 0) of these operators: 

'P = \Y^ i'jUx)^^,^^^, 0). (28) 

Here x, ip and 9 are spherical coordinates, j = |, |, 1^ = —j, — j + 1, • • • j, and z/ = e/x 
(where /i = ^J jU + 1) + A^ + | and e = ±1) are quantum numbers which label the eigenvalues: 

J = j(j + l)^>,!/5 ^2g^ 

Ql^j,K,U = ^^j,K,U- (30) 

The explicit form of ^j,K,u will be specified later. Notice that the eigenvalues of Qi can be 
found algebraically starting with relation (1211) and the first of relations (1291) . 

Substituting (12S1) into (1771) we obtain the following ordinary differential equations for radial 
functions ipj^u- 

-Q^ + -2 ] ^j., = E^j,,. (31) 



X^ X 
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This equation is solved by the following functions: 

^l^jnu = Cj^^x^-'^'^ exp ( ) ( a„ 2K + 1, 2\l -Ex ) (32) 



where Cj^v is an integration constant, J-" is the confluent hypergeometric function, and 



K=^p{p + l)-X\ a, = K + \-^=, u = ±\ 3{3 + l) + X^ + - (33) 



In order function fl32|) be bounded at infinity, the argument has to be a negative integer 



+ i ^ = -72, n = l,2,... (34) 

The corresponding eigenvalue ^ in (12 7p and the energy value E from ([2]) are given by the 
following equations: 

2 2 

a , „ ma 



where 



N = y/u{u+l) - + n+ n = 0,1,2,... (36) 

The energy levels (135|1 are degenerated w.r.t. quantum number k which are eigenvalues of 
the third component of the total orbital momentum. And there are no other degenerations. 
Let us present the basic functions ^j^K,v = ^j,K.u{f, which are used in formula 



if z/ = /i > 0, and 



if z/ = — /i < 0. Here ^j-i-i ^±1 are spherical functions. 
6.2 Supersymmetric system 

Let us discuss the maximally superintegrable system whose Hamiltonian H3 is defined by equa- 
tion f[T7|) . and consider the related eigenvalue problem: 

i/s^ = (-A + -cr ■ n + ^ - ^ = E^. (39) 

\ X X"^ X J 
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Equation (!39|) admits a number of symmetry operators given by relations ( !T3|) and ([Tj 
Among them there are two commuting integrals of motion J3, and two anticommuting 
constants of motion Qi and Q3. Symmetries J3 and make it possible to separate variables. 
The constants of motion Qi and Q3 enable to decouple the system of equations in radial 
variables. In addition, it will be shown that these constants of motion generate a specific 
degeneration of the energy spectrum. 

Like fl27|) . equation f l5I?]) is exactly solvable. Expanding its solutions via basis vectors fl?r|) 
and (1551) . i.e., using representation (1251) . we come to the following equations for radial functions: 

^-V',- = (-^ + ^^^^^ - f ) = M... (40) 

where u are parameters defined in (!33l) . 

Hamiltonian "H^, is shape invariant. In other words, it can be factorized: 

= a^a^ + Cy (41) 

where 

d d cP' 

= — + Wy, a+ = -— + Wy, Cy = -— —2, e = sign u, (42) 

dx dx (|i,| + i±i)^ 



and Wy is a superpotential: 

2a 2|z/|+e + l 



2|z/|+e + l 4x 



Moreover, operators iii, = Tiy — Cy satisfy the intertwining relations 

riyat = afHy+i. (43) 



This means that eigenvalues of Hamiltonian Tiy fl40|) can be found algebraically using tools of 
supersymmetric quantum mechanics. As a result we obtain these eigenvalues in the form fl35|l 
where 

N = n+\u\ + ^^, ra = 0,l,2, ... (44) 



Like in the Hydrogen atom, energy values fl35|) . fl44l) are proportional to the inverse square of 
the main quantum number N. But in contrast with the Hydrogen atom, N is not a linear com- 
bination of two independent non-negative integers, and so there is no an A'^-fold degeneration. 
However, since u can take positive and negative values as well, there is a two-fold degeneration 
typical for supersymmetric systems. In addition, since the quantum number k, does not affect 
the energy values, we have the additional (2j -|- l)-fold degeneration w.r.t. k. 

Let us show that the supersymmetric nature of spectrum (135|) . (14^ is caused integrals of 
motion Qi and Q3. 

It follows from (1211) that operator Q4 = ^[Qi,Q3] anticommutes with Qi and Q3. On the 
set of eigenfunctions of operators and Qi it is possible to define the rescaled operators 

Q3 = —-T^ V4 = — + , ^^ 45 
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which satisfy the following anticommutation relations 



QaQb + QbQa = 26abH. (46) 

^ 2 

Here any of subindices a and b independently takes the values 3 and 4, and H = + j2ffTy2- 

By construction both and (^4 commute with H, but this fact is also a consequence of f l46p . 

Thus the rescaled integrals of motion Qs, and Hamiltonian H form a basis of the su- 
peralgebra of supersymmetric quantum mechanics. The two-fold degeneration of spectrum of 
Hamiltonian H (and so of Hamiltonian i^a) is a direct consequence of algebraic relations (1461) . 

Notice that the ground state with energy E = -j^^ (see equations f l35|) . (jH]) for n = and 
e = —1) is not degenerated. Thus the supersymmetry of system (l39l) is exact. 

Let us also present the state vectors corresponding to eigenvalues (155]) . 



^P^^J^ = CnX^"^^'^ exp -J -Ex LF -n, 2|i/| + 1, a/ -Ex . (47) 



Like in (152]) . J^{—n, 2|z/| + l, x) is the confluent hypergeometric function. However, its arguments 
differ from arguments of function 



6.3 Equations for radial wave functions for arbitrary potentials 

Eigenvalue problem for Hamiltonians Hi and H2 can be effectively decoupled for the case of 
arbitrary function and / present in their definitions (1151) and f ITB]) . Here we deduce the 
decoupled equations for radial wave functions. 

For Hamiltonian Hi it can be done in complete analogy with section 6.1. Considering 
the eigenvalue problem for this Hamiltonian with arbitrary function (p and repeating all steps 
presented in this section before equation (jH]), we obtain the following equation: 

Hu^Jj^u = (—^+ ~ ^ + ifix)] = Ei)j^^ (48) 



where v is the parameter defined in 

Thus expanding solutions via eigenvectors of the commuting integrals of motion, satisfying 
fl29|) and fl30|) . it is possible to reduce the eigenvalue problem for Hamiltonian Hi f lTSj) . which is 
a three dimensional system of two coupled equations, to the infinite set of decoupled ordinary 
differential equations f HS]) . 

In order radial functions '4)jKu to have a good behavior at a; = 0, the potential ^{x) should 
increase not faster than ^ when x — )■ 0. If for small x this potential increases as then 
parameter a should satisfy the condition a > yl + A^ — |. 

For some particular potentials ip{x) equations ( HSj) can be solved explicitly. Examples of 
such potentials are: 

cfix) = --, (49) 

X 

y.(x) = - + 4, /3^0, (50) 
X x^ 

(p{x) = iJ^x^ . (51) 
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In the present paper only solutions corresponding to the Coulomb potential (H9l) are dis- 
cussed, see section 6.1. 

Consider now the eigenvalue problem for Hamiltonian H2 with arbitrary potential: 

W = (- V + o- ■ n/' + /2 - ^) ^ = Ei;. (52) 

In addition to (152|) . we impose on wave function ifj the following condition: 

= ((i^ ■ p + /) (o- • L + 1) + |o- ■ n) ^ = g7/> (53) 

where eigenvalues E and q are connected by the following algebraic relation: 

4g2 = (2j + lf E + c?. (54) 

Equations (1521) and (1531) are compatible since operators and Q2 commute each other and 
satisfy relations f|T6|) . 

To separate variables, let us expand solutions of equations fl52|) and fl53|) via the complete 
set of eigenfunctions of the commuting operators J^, J3, and cr ■ L + 1: 

Here ^^^^.^^((v', 6^) are spherical spinors, j, k, e = ±1 and /i = J + | are quantum numbers 
labeling eigenvalues of above mentioned operators: 

J ^j,K,epL J ( J ~l~ l)^j,K,e/^; 

(cr ■ L + 1) Qj^f^^i^f^ = efiflj^n,efj, 



As a result equation fl52|) is reduced to the following coupled system of ordinary differential 
equations for the radial wave function: 

+ M/i±l) +f-^) ,^^(a;) + /V,,.,-m(^) = ^^,,.,,(x), (56) 



dx"^ x"^ X 
d"^ /i(/i — 1) 2 



Using (jnSD this system can be decoupled. Substituting l^^ij into and turning into 
account the relations 

cr ■ p = cr ■ ncr ■ ncr ■p = (xp jcrnH — cr ■ n (cr ■ L + 1) , 



cr ■ n (cr • L + 1) = — (cr • L + 1) cr ■ n, 

1 

4' 

_ . d 

cr . nfi,- = fi,, X • pfi,- = 0, X ■ p^j,,^,,^ = -i^V',,.,./. 



fcr-L + D' = + i (o- -11)2 = 1, 
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we obtain the following system of first order equations: 

(/ - Q)'^j,K,,, + at'^j,K,-,, = 0, (58) 
-(/ + + = (59) 



were 



d fi a \ , f d fi a 



and 



'''' [dx^x 2fi)' [ dx^x 2/i 



^ 2j + l V (2j + l)2 

Solving equation ( 159|) for i'j^K-fi and substituting the derived expression into (l58i) or into 
(!56|) we obtain a decoupled system of second order equations: 

+ jq:^"^ + ~ j^j ^j''^''^ = ^^iA,/^- (60) 

Thus equation (13^ with arbitrary potential function / = /(x) can be reduced to the 
decoupled system of ordinary differential equations fl60|) for radial functions. However, it is 
doubtful whether the latter equations can be solved exactly for a fixed function / = f{x) if 
«//' ^ 0. 

To end this section, let us write system f l5B]) . flFTI) as a single equation with a matrix 
potential: 

= (-^ + k^,. = {e+ ^^A-M (61) 

where ^/ij^,^ = column(^/'j- V'j.k -m), and 

In other words, the effective potential U can be expressed via superpotential W . This circum- 
stance makes it possible to find the ground states ipfl ^ of Hamiltonians "H which are solutions 

of equation + W) V'j^i^^ = 0. The latter equation is nothing but the system (l58l) . ( 159|) with 
g = 0. 



7 Discussion 

We present the completed list of 3d Hamiltonians (|2]) with matrix potentials, which are invariant 
w.r.t. the rotation group and admit first and second order integrals of motion. This list is given 
by equations ( !T5|l - (fT8!l and includes four representatives, two of which are defined up to arbitrary 
functions. 
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The presented systems have a clear physical interpretation and describe particles with spin 
1/2 and non-trivial dipole moment. The system f llSp was discussed in paper [7]. Like the 
Hydrogen atom, it admits six integrals of motion satisfying algebra o(4). However, in contrast 
with symmetries of the Hydrogen atom, the Laplace- Runge-Lenz vector R in f|T8|) is dependent 
on spin. In the present paper we prove that this system is unique, and there are no other 2x2 
matrix potentials compatible with this symmetry. 

It is possible to verify by direct calculations that Hamiltonians H1..H4 do commute with 
the presented integrals of motion and do not admit another integrals of motion independent 
with ones given by relations (|T5|) - (|T8|) . A much more difficult problem is to prove that the list 
of superintegrable systems (fT5l) - (|T8l) is complete. To do this it is necessary to find all non- 
equivalent solutions of the determining equations (16|)- (|T2|) . The related detailed calculations 
can be found in the Appendix. 

Solving the determining equations we did not restrict ourselves to scalar and vector integrals 
of motion. In contrary, we consider also the tensor integrals of motion and prove, that for 
equation ([2]) they do not exist. 

The eigenvalue problems for Hamiltonians ([2]) are systems of coupled partial differential 
equations of second order, which are rather complicated. However, these Hamiltonians admit 
at least three commuting integrals of motion. Thank to this fact the eigenvalue problems can be 
effectively decoupled and reduced to ordinary differential equations in radial variables. Solutions 
for two of the eigenvalue problems which include the Coulomb potential are presented in section 
6, while the eigenvalues and eigenvectors of Hamiltonian (fT8|) have been found in paper [7]. 

In contrast with the case of scalar classical and quantum superintegrable systems, the sys- 
tems with spin have been classified only partially. Thus they belong to a perspective research 
field, and it would be interesting to extend our knowledge of these interesting and important 
subjects. In particular, it is desirable to classify second order superintegrable systems with 
dipole (and spin-orbit) interactions, which are not a priori rotationally invariant. This problem 
is much more complicated than in the case of scalar systems, and it is a challenge to solve it. 

One more interesting task is to search for relativistic counterparts of non-relativistic su- 
perintegrable systems. The relativistic analogues of the Pron'ko-Stroganov system [S] and the 
system with Hamiltonian ( 1T8|) are discussed in papers [2T] and [Ij, and it is not too difficult to 
extend this discussion to the cases presented in (|T5|) - (|T7|) . 

A Solution of the determining equation 

Here we present the main steps in solution of the determining equations (l6l)- (|T2l) . 
A.l Decoupling of the determining equations 

First we will show how the rather complicated system of determining equations (1^- (IT^ can be 
simplified and decoupled to algorithmically solvable subsystems. 

Let us start with the first order integrals of motion. The corresponding functions ^^"'^ are 
equal to zero, and equations ([6])- ([8]) are reduced to the following ones: 

^^a ^ ^ Q_ ^^^^ 
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These equations are easily integrated: 
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i 

where 

Here a^, v"", u, /i, fi'^, i/™" and /i™" are integration constants. Moreover, without loss of 
generality we can set a'' = since non-trivial correspond to already declared integrals of 
motion f|T3|l . 

Thus to find the first order integrals of motion it is sufficient to solve equations fl9|)- f|T2|) with 
the given coefficients flA2p and trivial ^i^"-^. Moreover, all cases enumerated in flA3p . flA4p and 
l5|) (which correspond to scalar, vector and tensors integrals of motion) should be considered 



separately. 

Consider the second order integrals of motion. In accordance with ([6j) functions should 
satisfy equations for Killing tensors of rank 2. Thus they are second order polynomials in 
which can be represented in the following form [20] 

where 

^O"'' = Ai^'^^ + X2{S^^x^ - x^x^), (scalar, even) (A7) 
•I.O"'' = X^^x'' + Xlx'' - 25"^A^Xc, (vector, odd) (A8) 

^Oab ^ ^ab ^ ^ab^2 _ ^ac^b^^ 
- X^'X^X^ + S'^'Xf 

^oa6 _ x^^c^^bd^^ ^ Xl^e^'^Xd. (tensor, odd) (AlO) 



(tensor, even) (A9) 



Here Ai, A2, A'^, A"^ are (real) integration constants. Moreover, A"^ are symmetric and traceless 
tensors. In the right brackets the covariant properties of integration constants and parities of 
as functions of indicated. 
Functions with a fixed value of m 7^ also are Killing tensors of rank 2 w.r.t. indices a 
and b. Their general form is analogous to (1A6P but more complicated thanks to the additional 
free index m: 

7 

^mab ^-^^mab (All) 
1=1 

where 

$r^ = A(2x'"r^-x"5'^^-x''5""), (scalar, odd) (A12) 
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+ (5™'^(x''A^x, - A^x^) + 5™^(x'^A^Xe - A^a;^) (vector, even) (A13) 

- x'"(25"^A^x, - Xlx'' - X^x^), 

^rnab ^ ^^mca^b ^ ^mcb^a^^^ ^ A^(5™"£^'='= + 5'^^e''''')Xc, (vector, odd) (Al4) 

if^mab ^mac\cb , ^mbc\ca , f ^rma ^dcb , iimb^dca\^ \dk^ 

(rank 2 tensor, even) (A15) 

(rank 2 tensor, odd) (A16) 



\ad m bdc \bd m adc 

J/ c -^c 4 X o J 



+ 2x™Af - {5'^''\\' + (5'"''A^")x„ 

^mab _ ^mac^cbd^^ ^ ^mbc^md^^^ ^^^^j^ 3 ^^^^^^^ Q^d) (A17) 

^mab _ ^ma6^2 _ ^mac^fe^^ _ ^mftc^a^^ ^ 5'^'' A^f '^'^x.Xd. (rank 3 tensor, even) (A18) 

Formulae (1A6P and flAllI) give the general solution of equations ([6]). One more subsystem 
of the determining equations which can be easy integrated is given by formula ([7j) . In this case 
we deal with the equation for Killing vectors, whose general solution is given by equation (lA2p . 

Solving of other equations (|^- (IT^ is a much more complicated problem which, however, can 
be effectively separated to relative simple subproblems starting with the following speculations. 

• By definition, Hamiltonian ([2]) admits integrals of motion (fT5|) and is an integral of motion 
by itself. Thus without loss of generality we can set = 0, $5"'' — ^'^d -^2^ (lA2p . 
flA6p and flA9l) respectively, since they correspond to higher order terms of operators if, 
Ja and JaJb- 

• Functions flA12p and flASp . flAlSp . flA14p correspond to scalar and vector integrals of 
motion while the remaining solutions generate tensor operators Since scalars, vec- 
tors and tensors transform in different way under rotation transformations which keep 
Hamiltonian invariant, all of them should satisfy the commutativity condition (jl]) inde- 
pendently. In other words, the determining equations (16|)- (|T2|) should be solved separately 
for scalar, vector and tensor operators. 

• Integrals of motion which are tensors of rank 3 are forbidden since the number of their 
components exceeds the maximal admissible number of integrals of motion. Thus it is 
possible a priori to set in flAllI) <l>™"'' = <l>^"'' = 0. 

• In accordance with f[T^ the external field F"- and potential are vector and scalar with 
well defined parities. Parities of solutions f lASp - flAlOp and flA12p - flA18p are transparent 
also. Then, analyzing properties of the remaining equations (18|)- (|T2|) under the space 
inversion we conclude that functions A™'" and Vt^ should have the same parity. In 
addition, the parity of functions (^^"■'^ ^ JSP°- and Vf* should be opposite to the parity of ^^"-^ ^ 
A™'" and QP^. Thus the system (jSD-dlSD should be solved separately for the cases when 
(^(^oab^ A°°, f2°) are even and odd with (<|)™°-''^ A™*^, fi™) being odd and even respectively. 

In accordance with the above the system of determining equations (18|)- (|T2|) is decoupled to five 
subsystems corresponding to scalar, vector and tensor functions (^^""^ with fixed parities. In 
other words we are supposed to solve these equations with given combinations of functions (^^^^ ^ 
^"^"^ and A°", namely: 

^Qab ^ ^mab ^ ^mab^ (Al9) 



14 



^Oab ^ ^Oab^ ^mab ^ ^mafe^ (A20) 

^Oab ^ ^mab ^ $mab^ ^Oa ^ ^ (A21) 

^Oafe ^ ^Oab^ ^mab ^ ^mab^ (A22) 

^Oa6 ^ ^Oab^ ^ma6 ^ ^maft (A23) 

where functions are defined by equations (lA8p - (lA10p and (lA12p - (lA16p . Except the case 
flA2ip the tensor A"'* should be trivial 

A. 2 Scalar constants of motion 

Consider first order scalar constants of motion. The corresponding coefficient functions A'^° are 
reduced to tensors A^" and A^" given in flA3p . The general form of scalar operator with 
such coefficient functions is: 

Q = z/io- ■ L + z/2cr ■ p + cr ■ ^f\x) + f{x) (A24) 

where and are arbitrary functions of x. The corresponding coefficient functions in ([5]) 
are: 

A'"'^ = z/i(5™" + z/2e"'"'=x„ ^° = /^ n"^ = x™/^ (A25) 

while functions ^'^"^ are trivial. 

Substituting (lA24p . (1A25P into equations (ITT]) and (IT^ we obtain that =Const and 
1^2 = 0. Then equation fllOp turns to identity, and the last remaining equation is solved by 
/i = ^ and fi = ^. As a result we obtain integrals of motion f fTSjl . 

Consider second order scalars. There exist the only second order scalar term specified by 
the coefficient function flA12p . The general form of the corresponding scalar operator is given 
by the following formula: 

Q = i(cr - pxL — cr-Lxp) + /(x)cr ■ L + (f{x)(T ■ p + g{x)(T ■ x + h{x) (A26) 

where f{x),(p{x), g{x) and h{x) are functions of x which will be specified in the following. 
The related multipliers A'"", Q"" and in ([5]) are: 

A'"'^ = (^(x)^'"'' + Q'^ = h{x), Q'^ = x"'g{x). (A27) 
Substituting f|T4|) . flA12l) and flA27l) into ([8]) we obtain the following conditions: 

^{x) = 0, fix) = x^i. (A28) 
The next equation, i.e., ([9]), is solved by the following functions: 

g{x) = -, F' = f{x)--. (A29) 

X X 

Equations (fTOj) and (fT2|) turn to identities, while the remaining equation ( ITTj) generates the 
following condition for function h{x): 

h\x) = 2xipi = 2f\x). (A30) 

Relations flA27p - flA30p together with ([2]), ^ specify infinite many Hamiltonians defined up 
to arbitrary function /. The corresponding integrals of motion are given by equations (|T6l) and 

IE 
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A. 3 Vector integrals of motion 
A. 3.1 First order integrals 

We start with first order vectors constants of motion. In accordance with (lA2p and flA4|) there 
are two possible sets of the corresponding non-trivial coefficient functions in ([5]): 

^ma^^mac^c^ l]° = Xe/i7l, = fi'' f2 + X^cX' fs (A31) 

and 

A°" = ai/^ A'"'^ = ^'^V^Xe - x'^z/^ ^^'^ = £'^'•^1/74 (A32) 

where /i,...,/4 are functions of x. Functions fi" and given in equations (lA3ip and (1A32|) 
represent generic scalars linear in vector parameters u and fi. Moreover, parities of these 
functions satisfy conditions requested in the last item in section Al. 

Substituting (lA3ip . (IA32p and (fn|) into determining equation we obtain the necessary- 
condition if = which correspond to the trivial external field F. Thus Hamiltonian with a 
nontrivial dipole interaction does not admit first order vector integrals of motion. 



A. 3. 2 General analysis of second order case 

Consider the second order vector constants of motion. The corresponding second order terms 
in ([5]) can include either even coefficients flA20p or odd coefficients flA2ip . 

Let us start with even functions f lA20p . Being included into the same integral of motion, 
all vector coefficients Aq, A4 should be proportional, i.e., = n^A" with some constants 
and A°. Thus we have the following coefficient functions for second order terms: 

$0'^'' = uoiX^x^ + AV - 26^^X^Xc), 

^"'^^ = uiX"'6'''' + z/2((5"'^A'' + S'^'^X'') + z/3(A™((5'*^x2 - x'^x^)) (A33) 
+ z/4(5""(x^A% - A^x') + 5'^\x''X'xc - A'x') - x"'{25''^X'xc - A'x'' - A^x")). 

The corresponding first order terms include even tensor functions A"^"" and odd vector func- 
tions A"'^. Moreover, all these functions should be linear in vector parameter A", and A°" should 
solve equation ([7]). Without loss of generality we can represent them in the following form: 

A°" = u^e^^^x^X", (A34) 

jyma ^ ^m^ack^c^kj^ ^ ^mka^k^c^cj^ ^ ^"^"^AVs (A35) 

where /i, /2 and /s are arbitrary functions of x. Moreover, it is possible to set 1/5 = 0, since a 
nontrivial z/5 corresponds to derivative terms of the total orbital momentum which is an integral 
of motion by construction. 

Vector function and scalar in ([5]) should be even and odd respectively. Their generic 
form is given by the following equation: 

= X^x^U = x^X^x^h + AVe- (A36) 

Considering the odd coefficients (lA2ip we have the following coefficient functions for second, 
first and zero order terms: 

= z/i(e'"'="A'' + e™'=^A'^)xc + V2X^ {S""" e^^^ + rt"'=^)xc, (A37) 
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A"*'^ = a;™x"A%7i + ^'""A^Va + x^A^s + XaX'^U, (A38) 
A"*^ = ^/3A^ = e^^^x^X^h (A39) 

where /1../5 are arbitrary functions of x. The other coefficients in ([5]) are trivial. 

Thus the classification of second order vector integrals of motion is reduced to finding non- 
equivalent solutions of equations (l6l)- (|T2l) for functions (]A33|) - (]A36[) and (]A37|) - (]A39p . 



mab 



A.3.3 Solutions for functions (K33l) -( IX36b with trivial $ 

For the case of trivial we have the following coefficient functions allowed by relations ([7]) 
and ([HD: 

while the generic form of fi" and $7° are still defined by equation flA36p . 
Substituting (fT4|) . (]A40[) and (1A36I) into equation ([9]) we obtain: 



( -(x'^A'' + x^'A'^) - 5''VA^ ) + - (xV'^a;'' - x°-x^X^rf) ^ 

Y 2 / 2 3!/ 



¥' = ::^' /5 = — , /6 = c (A42) 



-2z/(a;"A^ - 5"^A"x")<^ + AV^ + v'^x^h + S'^^x'X'f^ + x^x^X'^x'^ = 0. 

Equating coefficients for linearly independent terms x"A^, x^'A", x^x^'A'^x'^ and (5'^^A^x'^, we obtain 
the following system of equations: 

x^(l-4z/) + 2/^ = 0, (x^y + 2/5 = 0, 

y.' = 2/^, (2z/-l)y, + /5 = 0. ^ ^ 

It follows from f lAiTj) that 

1 a 2a 

^ = 7' = — ' /s = - 
4 x^ X 

where a and c are integration constants. 

The next determining equation, i.e., (fTOj) . reduces to the following form: 

(2z/F°' - ac)e'™"^xV'' = 0, 

or 

(A43) 

On the other hand, equation (fTTj) . which takes the form 
u'^xfr + xV^xV; + (z/'^x^ - xV=x")F°' = 0, 
is solved by the following functions: 

F' = h = - (A44) 

X 

where A is an integration constant. 

Relations f lA43p and flA44p are compatible if either a = OorA = c = 0. The first possibility 
corresponds to a trivial vector field F in equation If A = c = and a 7^ then formulae 
(]A42I) together with (]A40p . (]A36p . (|T4|) and ([5]) gives the Hamiltonian and integral of motion 
presented in equation (|T8|) . 
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A.3.4 Solutions for functions (lA33l) -( lA36b with nontrivial ^""''^ 

Functions (]A33p - (]A36p solve the determining equations (jSD and (j7]). The next step is to substi- 
tute flA33p - flA36l) and f|T^ into fl8l)- f|T2l) . Equating in ([8]) coefficients for linearly independent 
terms: 

ymab ^^mbk _j_ ^mak yb^^^k y^mab 2^^mck ^^c ^k ^ab y^mab 2,£^'^^ x'^ 

yrnab ^ ^^mbc^a ^ ^mac^b^^c^ ymafc _ ^rn^^a^bcd ^ ^b^acd^^c^d 

and using the identities 



(^^"'^bmk _|_ j^b ^amk\^^k yc^c ymab _|_ ymab ^2 ^mab _|_ ymab y^mafe-^ 



(A45) 



we obtain the following system of equations for functions /i, /s and /s: 

ViXip + - = 0, 

/2 - Z/3XV9 = 0, 

i^4xV + a;/i+a;V2 + /3 = 0, (A46) 

(z/gX^ - Vi)xip - X^f2 - /s = 0, 
Z/2XV9 - X(/2X)' + /3 = 0. 

The general solutions for system flA46p look as follows: 



^ ~ ( 2 ^ ^^ ' ~ / 2 , ^- ' 

/2 — — — — n"' /3 — TT; 



(A47) 



(x2 + a;)2 (x2 + u;)2 v^^ + v^ 

if z/3(z/3 + z/4) 7^ 0, z/2 = (z/4 + 2v^)uj, and 

V9 = a, /i = -az/i, /a = a(z/i + z/2), h = ^auix^ (A48) 

if 2/3 = 1/4 = 0. Here a is an integration constant. 

Consider solutions flA47p . The next determining equation, i.e., ([9]), is compatible with f lA47p 
with a nontrivial a only for the case a; = and z/i = z/2 = 1/3 = 0. In accordance with flA47p 
the corresponding Hamiltonian is reduced to the form given by equation f[T^ . Solving and 
the remaining equations (fTOj) and (ITTi) we obtain the following second order integral of motion: 

Q = J (cr ■ L + 1 + ao- ■ n) (A49) 

This result is trivial since (IA49P is a product of the known ffist order integrals of motion for 
Hamiltonian f lT5|) . 

Consider the case 2/3 = z/4 = and the corresponding solutions flA48p . Substituting functions 
flA33p . flA35p . flA36p and flA48p into and equating coefficients for linearly independent terms, 
we obtain the following system: 

2uoax - (2z/2 + z/i)x^ + z/gF"' + x/5 = 0, 

/g - z/Qttx + 2z/iF°' - z/ia^x^ = 0, ^^^^^ 
z/2F°' - Z/ott + x/5 = 0, 
2(z/i + z/2)a2x + = 0. 



This system has nontrivial solutions iff Xq = 0, i^i = —2^2, then 

h = -Ha'x' + 2ci), F° = + cix\ U = ^2 (^yx^ + Scix^^ + C2 (A51) 

where ci and C2 are integration constants. 

Relations flASip are compatible with equation f lTOj) only for V2<^ = 0. Thus a nontrivial 
second order symmetry (which corresponds to z/2 7^ 0) is possible for the trivial field F. We 
show that, except flA49p . there are now second order integrals of motion with even functions 
(]A20|) provided ^'""^ are nontrivial. 



A.3.5 Solutions for functions (JMB-dMHD 



Functions ( ]A37|) -( lA39|) solve the subsystem of the determining equations given by relations 



and ([7]). Substituting these functions into the next equation ([8]) and equating the coefficients for 
linearly independent terms x"'x''x"^X'^x'^, x'^S"'''x'^\'^, {5'^"'x'' + S'^^x"')x'^X'^, X'^x'^x^, x'^{X°'x'' + 
A^x'^), X'^5^'' and 5""'X'' + 5'^^X'' we obtain the following system: 



(A52) 



/i = T^2V, fi = 0, 

f2 + 2/ix = 0, + (/i + (z/2 - ^ = 

fi + I/2XV = 0, /4 - V2X^ = 0, 
/2 + /s + (Z^l - Z^2)xV = 0. 

Its general solutions are 

fl = 14 = 1^2= 0, /2 = /i, /s = /i, f= ^ (A53) 

where a and /i are integration constants. 

The next determining equation takes the following form: 

^c^c^abk^k + AVs + x^e^^^X^x'^ = 0. 



Its solutions are = 0, F° = — ^ + Const. However, these solutions and solutions f lASSp are 



compatible with equation f[T2l) only for a = 0, i.e., for the case of the trivial Pauli interaction. 

Thus we prove that vector constants of motion exist only in two cases, namely, when the 
Hamiltonian is given by equation f[T^ or f[TB|) . In the first case integrals of motion are given by 
equation flA49p and so are products of the first order symmetry operators. In the second case 
the integrals of motion are components of the Laplace- Runge-Lenz vector with spin, introduced 
in [3. 



A. 4 Tensor integrals of motion 

In accordance with the analysis presented in section Al, the tensors of the zero and first order 
in do not exist. Moreover, to search for second order tensor integrals of motion we can a 
priori restrict ourselves to the case when the coefficients for second order terms are given either 
by equation flA22p or by equation flA23p . Consider both versions consequently. 
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A. 4.1 Integrals of motion with odd and even <|)™'^'' 

The corresponding coefficient functions for the second order terms are given by equations (lA22p , 
flAlOl) and flA15l) . where the second order constant tensors should be proportional: 

$0"'' = uiiX'^'e''''' + A^%"'='^)x^ (A54) 

^mab ,, /^mac\cb , ^mbc\ca\ 

+ z/3 + rt"'=^)A'=Vx'^ - a;"(A''^e^*^" + A^*^e''^")x") . 

The related tensors A™" and vectors Q"^ should be even functions of x linear in tensor 
parameter A'*''. All such functions can be represented in the following form: 

^ma ^ ynaj^ ^ x'^A'^X^a + x'^X^'x'^h + 5""'A'='^X (A56) 
_ ^mkc^cb^k^bj^ ^^57) 

where /1../5 are arbitrary functions of x. 

The coefficient should be a scalar odd function linear in A"''; such functions do not exist. 

Substituting flASSp flA56p and f[T^ into ([H]) and equating coefficients for the linearly indepen- 
dent terms x'"(5"^A'=Vx'^, A'^'^x'' + A'"^x^ (5'"''A^'=x^ + 5'^^A'^'=x'=, x'"A'^^ ^'^^x^ + ^'^V, (5"^A'"'=x'=, 
x"x''A™''^Xc and x'^(x^A"'^ + x^A*"^), we obtain the following conditions: 



(A58) 



1^3 = 0, /2 = 0, /s =0, /4 = /i, 

+ U2Xip = 0, 2/4 - Z/2V3 = 

where is an integration constant. 

Let 7^ then it follows from flASSp that 

2/yr™' 

= — , A™'' = /i(5'"*A"Vx'^ - A'"''x^). (A59) 

If 1/2 = then /4 = 0, /i = /i, and so 

$""^^ = 0, A'^^ = /iA™^ F"^ = x™(^ (A60) 

where /i is a constant and is an arbitrary function. 

Substituting flA54l) . flA59p and flA55l) with 1/3 = into the next determining equation ([9]) 



and equating to zero the coefficient for the term ^'^"■k^k^cd^c^d obtain the condition /i = 0. 
In accordance with flA59p it corresponds to the trivial external vector field and so to the trivial 



Pauli term in Hamiltonian ([2]). 

Consider now the version represented by equations flA60p . Substituting flA54p and flA60p 
into dS]), using the identities 



~ ' fA61) 

^mak^k^cd^c^d ^ r^rn^akc^kycb^b _ ^mkc ^k ^cb ^b _ ^^^macyck^k 

and equating the coefficients for linearly independent terms presented in the r.h.s. of equation 
(lA6ip we obtain the following conditions: 

</.' = 0, /^ = 0, 2z/iyp-/5 = 0, i/i</p - 2/iyp + /s = 0. (A62) 
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It follows from ( 1M2|) that 



^ = h = 2z/z^i, /i = -ui (A63) 



where z/ is a constant. However, equations flA57l) . flA60p and flA63l) are incompatible with ([10 



ii (f = v ^ 0, thus there are no tensor integrals of motion of the considered type for Hamiltonian 
([2]) with a non-trivial matrix term cr ■ F. 

A. 4. 2 Integrals of motion with even and odd 

The related coefficient functions for the second order terms are given by equations flA23p . flA9p 
and (]A16I) . Moreover, without loss of generality it is possible to set in f lAOp = since the 
terms including this coefficient can be related to the obvious symmetry operators Xf'JaJb where 
Ja and Jb are components of the total orbital momentum. Thus we have: 

$0"'' = jjiX^b (A64) 
= z/2(A"^'^a;^ + A'^V - 25'^''A'""xc) + U3{2x"'X''^ - (5'^'^A^" + 5'"^A'^")xc). (A65) 

The corresponding tensors A™" and vectors fi™ should be odd functions of x linear in tensor 
parameter A*^^. The generic form of such functions is given by the following equations: 

^mb ^ ^rnchyck^kj^ ^ £^^'=A'"V2 + e"'''' X^^ x'^ h + ^'"''"A^ V/4, (A66) 

_ x^'^x'^h + x^A'^Vx'^/e (A67) 

where fi..f& are functions of x. 

In addition, there is a possible odd non- derivative term fi" in ([5]): 

Substituting (1A65|) and (1A66P into ([8j) we obtain the following equation: 



(A69) 



Xmab f 2 , \rmab f 1 i \^mab f 1 , r) -irmab r , \rmab f 3 i r 
X X^ X X 

f 

( Vtnab I vmab\ r , vmab J 4 /,, \rmab ,, \rmab , ,, vmab\ „ n 
- (Ai + A2 )/4 + — - (Z/2A7 - Z/aAg + Z/3A5 )V5 = 

where 

j^rnab ^ ^mk^^a^bck ^ ^b^ack^^c^ j^mab ^ ^m^^ak^bck ^ X^k^ack-^^c^ 

j^mab ^ ^^ma^bdk ^ ^rnb^adk^^d ^kc^c ^ j^mab ^ ^S^^ e"''''')x^ XV , (A70) 

j^mab ^^mbcyak _j_ ^^'^'^ X^^^x'^X^ j^mab ^^mcbyac _j_ ^mca ybc^^^2 

j^mab ^ ^^a^mbc ^ rjfi ^n^ac^X^d ^k ^d ^ ^mab _ 2x^ X^nickx" X^'^ x'^ . 

Tensors (1A70P satisfy the following conditions: 
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-mab 
3 


j^mab 


vmab 
- ^3 




■mab 
4 


j^mab 


1 v"mab 
+ -^8 


j^mab 


■mab 
5 


j^mab 


1 vmab 
+ ^6 


1 j^mab 



while all terms in the r.h.s. of (IA7ip are linearly independent. Thus equation (1A69P generates 
the following system of equations: 

/i = 0, /^ = /I = o, 

2/3 = {1^2 + y^)^, h = J^sV, (A72) 

/2 = UX, /2 + fiX + UsipX = 0. 

Solutions of this system with a non-trivial cp are: 

fi = U = 1^3 = 0, /2 = /i, /3 = ^, ^ = \ (A73) 

where fi and A are arbitrary constants. 

In accordance with (]A73p coefficient functions (1A66P and vector field F° are reduced to the 
following form: 

Substituting (IM4D . (IMSj) with z/3 = 0, (IMTD and flX74D into © and equating the coeffi- 
cients for linearly independent terms X°'^,6°''^X'^'^x'^x'^,x°'x"^,x'^X"'^x^ and x°'X^^x^ we obtain 
the following system: 

= 0, 
z/iA + /5 = 0, 

/e + ^^^A^x^ - = 0, 

/e - ^2A=^x = 0, 
z/2F°' + + 2/63; = 0. 

This system is compatible only for z/2 = or A = 0. The first case is compatible with flTU]) 
only for A = 0, thus all admissible solutions correspond to a trivial external vector field F"-. 

Summarizing, there are no rank 2 tensor integrals of motion for equation (|2]) if F 7^ 0. 

Integrals of motion which are tensors of rank 3 can be excluded a priory since they have 
too many independent components whose number exceeds the maximal possible number of 
constants of motion for a 3d system. The absence of such integrals of motion (and integrals 
being tensors of rank R>3) also can be proven directly using the determining equations (P)- (IT^ . 
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